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During  the  three  and  half  years  of  the  tenure  of  this  grant,  we  constructed  our 
research  program  and  carried  out  a  number  of  important  studies,  as  summarized  below.  In 
particular,  we  wrote  microcanonical  (constant  NVE)  and  canonical  (constant  NVT) 
ensemble  MD  codes  for  studying  processes  on  the  (100),  (111),  or  (1 10)  faces  of  any  fee, 
bcc,  or  diamond  lattices.  We  also  wrote  an  isobaric-isothermal  (NPT)  ensemble  MC 
program  to  study  equilibrium  properties  of  SiGe  alloys  and  superlattices.  Although  true 
deterministic  dynamics  is  only  available  in  the  NVE  ensemble,  most  surface  science 
experiments  are  carried  out  with  conditions  of  constant  temperature,  constant  volume,  and 
constant  number  of  particles  (NVT).  We  chose  to  operate  in  a  mixed  NVE/NVT  ensemble, 
where  the  near  surface  layers  were  treated  with  correct  dynamics  (NVE)  and  the  bottom 
layers  were  treated  with  the  canonical  ensemble  (NVT)  in  order  to  mimic  the  actual 
isothermal  reaction  conditions.  The  NPT  ensemble  is  important  to  utilize  in  cases  where 
the  lattice  constants  of  the  alloy  are  unknown  (i.e.,  the  volume  of  the  crystal  is  allowed  to 
fluctuate). 

As  a  first  step  toward  studying  adsorbate  dynamics,  we  studied  the  reconstruction 
dynamics,  equilibrium  structure,  and  short-time  dynamical  behavior  of  the  bare  Si(100)  and 
Ge(100)  surfaces.1’2  In  particular,  we  carried  out  finite  temperature  MD  simulations  using 
empirical  two-  and  three-body  interaction  potentials  for  Si  and  Ge  that  have  been  fit  to  bulk 
properties.3  Our  simulations  reconciled  several  seemingly  conflicting,  temperature- 
dependent  surface  science  experiments  in  terms  of  a  dynamical  model  where  buckled  and 
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symmetric  surface  dimers  interconvert  on  a  subpicosecond  time  scale  via  a  surface  phonon 
(the  dimer  buckling  vibration).  In  particular,  the  structure  of  the  surface  observed  was 
shown  to  depend  on  the  time  scale  of  the  probe.1-2  We  were  able  to  explain  both  long  time 
scale  (scanning  tunneling  microscopy)  and  short  time  scale  (electron  diffraction  and 
photoemission)  experiments  based  on  this  dynamical  model.  On  a  technical  note,  our 
mixed  NVT/NVE  ensemble  was  found  to  produce  much  more  ordered  surface  structures 
than  NVE  dynamics  did.2  Thus,  it  is  important  to  use  this  mixed  ensemble  during  all 
dynamical  simulations  that  involve  evolution  or  absorption  of  heat  (NVT  required)  and 
where  time-dependent  information  is  desired  (NVE  required). 

We  performed  isothermal-isobaric  ensemble  (NPT)  MC  calculations  of  the 
structures  of  Si-Ge  alloys4  and  superlattices5  as  a  function  of  composition  and  temperature. 
We  predicted  the  equilibrium  structures  of  the  bulk  alloy,  the  alloy  surface,  and  Si-Ge/Si 
and  Si-Ge/Ge  superlattices.  (Note  that  the  superlattice  work  is  still  in  progress  and  is 
expected  to  be  completed  this  summer.)  Since  diffusion  of  atoms  in  a  bulk  lattice  is 
extremely  difficult,  we  developed  an  algorithm  which  speeds  up  the  equilibration  process 
considerably.  We  carry  out  conventional  MC  random  walks,  but  every  500  steps  we 
attempt  to  permute  the  positions  of  two  atoms  picked  at  random.  Then  we  allow  the  nearest 
neighbors  to  relax  around  atoms  that  may  have  switched  and  then  continue  MC  steps  for 
another  500  moves.  We  find  that  the  main  difference  between  the  bulk  structure  (which  we 
found  to  be  a  random  alloy  above  =150  K)  and  the  surface  of  an  Si-Ge  alloy  is  that,  given 
the  opportunity,  Ge  will  segregate  to  the  surface.  This  occurs  because  the  Ge  surface 


energy  is  lower  than  that  of  Si  (i.e.,  Ge-Ge  bonds  are  weaker).  This  has  important 
implications  for  epitaxial  growth  of  such  alloys,  in  that  if  surface  segregation  is  to  be 
avoided,  growth  should  be  carried  out  at  a  low  enough  temperature  to  kinetically  limit  this 
surface  segregation.  Features  of  the  superlattice  simulations  include  identification  of  a 


mechanism  for  relieving  strain  at  the  interface  between  the  alloy  and  the  pure  element:  the 
alloy  expands  or  compresses  in  the  direction  normal  to  the  interface.  These  alloy 
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simulations  are  the  first  step  toward  studying  the  structures  of  analogous  systems  such  as 
SiNx/Si  and  SiFx/Si  interfaces. 

Our  first  forays  into  studying  silicon  surface  chemistry  involved  modelling 
chemisorption  of  H  and  F  on  Si(100)-2xl  using  first  principles  quantum  mechanical 
calculations  on  embedded  silicon  clusters  (SixHy,  x=2-9,  y  =  2-12)  chosen  to  represent 
truncated  pieces  of  the  Si(100)  surface.6*9  The  embedding  procedure,  which  surrounds  the 
subsurface  Si  atoms  with  a  bulk-like  environment  is  described  elsewhere.8  We  have 
developed  several  new  theoretical  strategies  for  obtaining  reliable  adsorption  energetics:  we 
use  large  clusters  to  obtain  analytic  gradient-optimized  geometries  at  the  generalized  valence 
bond  (GVB-PP)  level10  and  then  map  those  geometries  onto  smaller  clusters  (e.g.,  two  to 
six  silicon  atoms)  when  performing  correlation-consistent  configuration  interaction  (CCCI) 
calculations  to  obtain  bond  energies  and  heats  of  reaction.11  Based  on  the  largest  Cl 
calculations  we  can  do  on  the  very  large  clusters  (essentially  a  complete  active  space 
(CAS)-CI),  we  estimated  the  error  incurred  by  using  the  smaller  clusters  to  represent  the 
surface  to  be  =0.1  eV  8  In  order  to  properly  model  chemisorption,  it  is  necessary  and 
sufficient  to  describe  the  local  wavefunction  accurately;  we  have  shown8  that  choosing  the 
optimal  geometry  for  the  smaller  cluster  to  be  that  for  the  larger  cluster  leads  qualitatively  to 
the  same  local  surface  wavefunction  for  both  clusters. 

Since  the  cluster  model  of  chemisorption  normally  treats  the  limit  of  zero  coverage, 
which  is  clearly  unrealistic  for  most  situations,  we  also  developed  a  means  of  estimating 
binding  energies  for  chemisorption  at  high  coverages.8  Essentially,  this  involves 
evaluating  the  lateral  interactions  between  coadsorbates.  In  particular,  we  calculate  the 
interaction  energy  (normally  nonbonded  repulsions)  between  adsorbates  on  two  small 
clusters  as  a  function  of  distance  between  the  clusters.  The  atoms  in  these  clusters  are 
oriented  in  such  a  way  so  that  the  dominant  contributions  to  nonbonded  repulsions  are  due 
to  the  adsorbates.  These  repulsive  forces  are  evaluated  at  the  GVB-PP  level,  which  should 
provide  an  upper  bound  on  the  repulsion  energy  (since  dynamical  correlations  not  included 
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in  this  wave  function  would  tend  to  decrease  this  repulsion).  We  are  careful  to  evaluate 
counterpoise  corrections  as  a  function  of  distance  for  basis  set  superposition  errors  that  are 
often  present  in  weakly  interacting  systems.  Since  we  calculate  this  repulsive  interaction  as 
a  function  of  distance,  we  are  able  to  use  a  criterion  of  minimal  repulsions  to  establish 
equilibrium  adsorbate-adsorbate  distances  on  the  surface.8 

The  above  techniques  were  used  to  predict  heats  of  adsorption  for  H  atoms  in  the 
monohydride  phase  on  Si(100).6  Understanding  hydrogen  interactions  with  Si(100)  is  a 
critical  first  step  before  undertaking  studies  of  nitridation  of  silicon  by  NH3  or  oxidation  of 
silicon  by  H2O,  since  H  atoms  are  generated  in  both  of  those  reactions  and  desorption  of 
hydrogen  is  in  fact  a  rate-limiting  process.  We  predicted  that  the  Si-H  bond  strength  is 
=3.8  eV  and  that  the  endothermicity  for  desorption  of  H2  from  this  phase  is  =3  eV,  which 
is  significantly  higher  than  any  of  the  desorption  activation  energies  measured  by  LITD  or 
TPD.12d3  Only  two  ways  exist  for  the  measured  activation  energy  to  be  lower  than  the 
endothermicity:  either  the  measured  activation  energy  is  an  "apparent"  activation  energy, 
which  is  the  difference  between  (at  least)  two  activation  energies  for  competing  steps,  or 
the  pathway  for  hydrogen  desorption  does  not  involve  the  species  for  which  the 
endothermicity  was  calculated.  We  first  suggested  that  desorption  and  diffusion  might  be 
those  competing  steps.6  Thus,  we  also  calculated  activation  energies  for  surface  diffusion 
of  H  on  Si(100)-2xl.9  Several  different  pathways  for  surface  diffusion  were  studied:  (i) 
hopping  from  one  dangling  bond  (db)  to  another  on  the  same  Si  dimer,  (ii)  hopping  from 
one  db  to  a  db  on  an  adjacent  Si  dimer  in  the  same  row;  and  (iii)  hopping  from  one  db  to  a 
db  on  an  adjacent  Si  dimer  in  the  adjacent  dimer  row.  Our  multiconfiguration  self- 
consistent-field  and  Cl  calculations  for  these  diffusion  pathways  predict  that  the  barrier  is 
lowest  for  hopping  from  one  db  to  a  neighboring  db  on  an  adjacent  dimer  in  the  same  dimer 
row.  This  makes  sense  upon  reflection.  For  example,  hopping  from  one  db  to  another  on 
the  same  dimer  involves  moving  across  a  bridging  transition  state  that  has  small  overlap 
with  either  db,  since  both  db’s  point  away  from  the  bridging  position.  Hopping  from  one 
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db  to  a  db  on  an  adjacent  Si  dimer  in  the  next  row  over  is  bad  because  the  Si  atoms  between 
which  the  H  atom  is  hopping  are  = 5.3  A  apart;  thus,  the  orbital  overlap  is  low  again  in  the 
transition  state.  The  reason  the  barrier  is  lowest  for  hopping  between  dimers  in  the  same 
row  (in  the  [lIO]  direction)  is  because  the  orbitals  between  which  the  H  atom  is  hopping 
are  only  3.84  A  apart .  Thus,  there  is  a  preferred  direction  of  diffusion  along  dimer  rows! 
Quantitatively,  the  barrier  for  this  most  favorable  pathway  is  predicted  to  be  about  2.1  eV. 

Recently,  we  have  carried  out  calculations  of  various  possible  H2  desorption 
pathways  for  both  the  monohydride  H-Si-Si-H  and  dihydride  H-Si-H  species  on  Si(100), 
as  a  function  of  coverage.14  We  find  that  the  barrier  (calculated  at  the  CASSCF  level)  for 
direct  desorption  from  the  monohydride  state  (a  1,2-elimination  reaction)  is  over  100 
kcal/mol  and  hence  is  ruled  out  as  a  mechanism  for  desorption  of  H2  for  submonolayer  or 
monolayer  coverages  of  hydrogen  on  Si(100).  This  originally  led  us  to  propose  as  a 
mechanism  hydrogen  desorption  from  two  adjacent  monohydrides  in  the  same  dimer  row,6 
but  this  direct  mechanism  still  suffers  from  the  fact  that  it  will  lead  to  a  barrier  at  least  as 
high  as  the  reaction  endothermicity,  which  we  calculated  to  be6  higher  than  the  observed 
barrier.  Thus,  as  mentioned  above,  a  direct  mechanism  is  out  of  the  question.  We  had 
suggested,  as  mentioned  above,  that  the  barrier  might  be  the  difference  of  two  activation 
energies  of  competing  steps,  desorption  and  diffusion,  but  further  analysis  of  the  kinetics 
of  such  a  mechanism  reveals  that  this  is  not  feasible.  Thus  we  were  forced  to  conclude  that 
desorption  at  monolayer  and  submonolayer  coverages  of  hydrogen  on  Si(100)  could  not  be 
due  to  desorption  from  either  two  H-Si-Si  or  one  H-Si-Si-H  groups.  The  answer  came 
when  we  examined  the  desorption  activation  energy  for  the  dihydride  species,  SiH2,  and 
found  that  it  varied  from  60  kcal/mol  at  low  coverage  to  43  kcal/mol  at  2  ML  coverage  of 
hydrogen.  As  it  turns  out,  the  latter  value  is  in  good  agreement  with  thermal  desorption 
data  for  the  dihydride  desorption  state  at  high  coverage.13  The  dihydride  state  at  2  ML 
coverage  of  hydrogen  desorbs  at  a  temperature  w  ell  below  that  of  the  monohydride  state, 
presumably  because  the  higher  surface  coverage  induces  lateral  repulsions  (we  calculate 
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them  to  be  17  kcal/mol  at  2  ML  coverage)  that  destabilize  SiH2  on  the  surface. 
Interestingly,  the  low  coverage  limit  for  the  activation  energy  for  H2  desorption  from  SiH2 
is  in  good  agreement  with  the  activation  energy  observed  for  desorption  of  H2  at  low 
coverages,  which  was  thought  to  proceed  via  the  monohydride  state.  We  are  now 
convinced  that  the  desorption  of  hydrogen  at  all  coverages  up  to  2  ML  proceeds  via  SiH2 
and  that  the  difference  in  activation  energies  between  di-  and  mono-hydride  states  is  merely 
due  to  differential  lateral  repulsions  at  different  coverages.  This  is  consistent  with  the 
observed  first  order  kinetics,  the  measured  activation  energies,  as  well  as  with  laser- 
induced  fluorescence  studies  by  Zare  and  co-workers15  which  concluded  simultaneously 
with  our  work  that  desorption  of  H2  from  the  mono-  and  di-hydride  states  goes  through  a 
common  intermediate,  which  they  suggest  to  be  SiH2.  In  their  work,  they  suggested  a 
possible  isomerization  mechanism  whereby  H-Si-Si-H  converts  to  H2Si  +  Si,  in  order  to 
produce  the  SiFb  intermediate.  Our  original  proposal  was  that  H  diffusion  to  unsaturated 
HSi  defect  sites  might  be  the  mechanism  for  formation  of  Sifb  at  low  coverages  (since  we 
calculate  the  diffusion  barrier  to  be  much  lower  than  any  desorption  barrier,  it  will  not  be 
rate-limiting).  This  could  still  be  happening,  but  we  have  decided  to  calculate  the  barrier 
for  their  mechanism  to  see  if  it  is  also  energetically  feasible;  this  wwk  should  be  completed 
this  month  and  then  we  will  be  submitting  a  u'hole  story  on  hydrogen  desorption 
pathways.14 

Turning  to  another  area,  w'e  have  used  the  above  quantum  mechanical  methods  to 
predict  adsorbate  structures  and  energetics  of  the  initial  stages  of  F-atom  etching  (i.e.,  F 
adsorption  and  reaction)  of  Si(100)-2xl  at  0  K.7>8  We  find  that  one  monolayer  (1  ML)  of 
F  atoms  adds  easily  to  this  surface  without  breaking  any  Si-Si  bonds,  releasing  a  large 
amount  of  energy  (=6  eV/Si-F  bond  formed).  Above  1  ML,  it  is  exothermic  to  form  a 
surface  that  is  75%  SiF(a)  and  25  %  SiF2(a),  but  then  2.5  eV  endothermic,  if  carried  out 
sequentially,  to  form  a  surface  that  is  half  SiF(a)  and  half  SiF2(a).  Above  this  coverage  (1.5 
ML),  we  find  it  is  very  difficult  to  add  more  fluorine,  because  of  large  repulsions  between 
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the  F  atoms.  This  is  consistent  with  both  molecular  beam  studies,  which  found  that  the  rate 
of  fluorine  adsorption  decreases  dramatically  above  1.5  ML,16  and  XPS  studies  that  find  a 
stable  phase  is  formed  that  is  close  to  3:1  SiF(a)  to  SiF2(a).17  Our  calculations  are  the  first 
to  provide  such  thermodynamic  data  and  to  predict  a  specific  fluorinated  surface  structure. 

We  then  used  our  ab  ;nitio  CCCI  data  for  F  interacting  with  Si  and  other  F  atoms  to 
refit  the  parameters  of  an  analytic  potential  function  for  Si-F  due  to  Stillinger  and  Weber 
(SW).18  First,  we  compared  predictions  of  the  SW  empirical  potential  (which  was  fit  to 
the  only  Si-F  data  available,  that  of  the  gas  phase  SiFi .4  series)  to  all  of  our  ab  initio 
reaction  paths.19-20  We  showed  that  the  SW  potential  yielJs  very  large  errors  (up  to  1-4 
eV  for  binding  energies)  when  compared  to  our  cluster  results.  The  SW  function  is  simply 
too  repulsive  for  nonbonded  interactions  and  is  not  attractive  enough  for  bonding 
interactions.  Indeed,  this  explains  why  the  two  sets  of  constant  temperature  MD 
simulations  using  this  potential  found  only  saturation  of  the  Si  dimer  db's;  no  etching  or 
breaking  of  dimer  bonds  was  observed.21*22  Clearly,  the  SW  empirical  potential  is 
insufficient,  since  it  does  not  predict  that  fluorine  etches  silicon  at  room  temperature,  when 
in  fact  it  does  so  spontaneously.  However,  since  the  overall  morphology  seemed  correct, 
we  retained  almost  the  same  SW  functional  form  and  searched  for  new  parameters  to  fit  our 
ab  initio  data.  We  developed  a  simulated  annealing  algorithm  to  sample  the  12-parameter 
space  of  the  Si-F  two-  and  three-body  terms  in  an  unbiased  manner,  using  =50  pieces  of  ab 
initio  data.19-20  We  then  programmed  the  forces  for  our  new  potential  designed  to  properly 
describe  the  fluorosilyl  layer.  We  then  performed  the  first  ab  initio-derived  dynamics  for  a 
surface  chemical  reaction,  selected  results  of  which  we  now'  describe.19-20 

We  exposed  an  equilibrated  Si(  1  ()0)-2x  1  surface  to  a  fairly  high  pressure  of  F 
atoms  (one  F  atom  was  added  a  few'  A  above  the  surface  every  62.5  fs,  which  corresponds 
to  p  =  200  atm)  up  to  a  coverage  of  2.3  ML  of  fluorine.  Two  sets  of  NVT/NVE  trajectories 
with  identical  initial  conditions  (Tsurfacc  =  300  K  and  the  F  atoms  had  thermal  velocities 
corresponding  also  to  300  K)  were  followed,  one  as  a  control.  We  attempted  to  etch  Si 
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with  our  potential  in  one  trajectory  and  with  the  original  SW  potential  in  the  control 
trajectory.  We  did  this  because  the  reported  MD  simulations  with  the  SW  potential  that  did 
not  etch  Si  were  done  at  =20  atm  pressure18  and  we  wanted  to  make  sure  that  any 
differences  in  the  two  trajectories  were  rigorously  due  to  the  change  in  potential  and  not  the 
change  in  pressure.  Indeed,  we  find  that  the  SW  potential  still  does  not  etch  Si  even  at  this 
pressure  that  is  lOx  higher  and  still  only  saturates  the  db's  on  the  Si  dimers.  By  contrast, 
we  find  that  our  potential  reacts  readily  with  Si  at  room  temperature,  producing  SiF,  SiF2, 
and  SiF3  on  the  surface.  However,  contrary  to  the  0  K  quantum  calculations  described 
previously,  the  heat  generated  by  the  formation  of  6  eV  Si-F  bonds  leads  rapidly  to  a 
disordered  phase.  Further,  fluorine  atoms  attack  not  only  the  Si  atoms  in  the  surface 
dimers,  but  also  attack  subsurface  atoms.  At  low  coverages,  we  find  that  the  F  atoms 
mostly  adsorb  onto  the  Si  surface  dimer  atoms  along  the  db  direction  (=23°  off  the  surface 
normal)  and  hence  are  merely  saturating  the  db’s.  However,  even  at  0.5  ML  of  F  atoms, 
we  find  that  SiF2  forms  on  the  surface  as  a  minority  species.  As  more  and  more  F  atoms 
adsorb  onto  the  surface,  we  find  that  the  near  surface  layers  become  disordered,  and  we 
further  find,  based  on  our  adsorption  probability  studies  discussed  below,  that  this  disorder 
is  key  for  etching. 

No  mechanistic  models  of  silicon  etching  have  been  developed,  other  than 
macroscopic  kinetic  models  that  lack  microscopic  details.23  Our  simulations  have  already 
produced  information  that  was  unavailable  from  experiments,  e.g.,  that  fluorine  atoms 
attack  not  just  the  dimers  at  1  ML  coverage  but  also  subsurface  Si  atoms  and  that  even  at 
this  coverage,  the  huge  heat  release  per  Si-F  bond  formed  leads  rapidly  to  a  disordered 
surface  and  near-surface  structure  with  many  broken  Si-Si  bonds.19  As  mentioned  above, 
these  trajectories  also  revealed  that  induced  disorder  may  be  key  to  the  etch  rate.  We 
pursued  these  ideas  by  comparing  the  reaction  probabilities  for  perfectly  ordered  overlayers 
of  adsorbed  fluorine  to  those  for  disordered  surfaces.20  We  find  that  while  the  adsorption 
probability  for  atomic  fluorine  is  unity  on  the  clean  Si(l()0)  surface  at  room  temperature,  it 
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is  zero  on  a  perfectly  ordered  1  ML  surface  of  fluorine  (i.e.,  where  the  surface  is  covered 
with  F-Si-Si-F  groups).  If  we  instead  consider  a  disordered  1  ML  surface  that  is  produced 
via  the  above  etching  simulations,  we  find  that  the  adsorption  probability  is  significant 
(0.15).  Thus,  it  appears  that  ordered  overlayers  (e.g.,  repeating  units  of  F-Si-Si-F  dimers) 
inhibit  etching  dramatically.  In  other  words,  the  reaction  is  dependent  on  the  continual 
production  of  more  and  more  dangling  bonds  for  F  atoms  to  attack,  which  is  the  case  for 
the  disordered  surface.  Further,  we  find  that  ordered  surface  structures  above  1 .25  ML 
are  unstable  with  respect  to  highly  disordered  surfaces.  Thus,  we  have  concluded  that  this 
reaction  proceeds  in  particular  by  local  heating-induced  disorder,  where  the  local  heating  is 
due  to  the  highly  exothermic  formation  of  Si-F  bonds  and  where  disorder  catalyzes 
formation  of  dangling  bonds  necessary  for  the  etching  to  proceed  onward. 

As  it  is  always  the  case  that  one  wants  to  treat  more  and  more  realistic  and  hence 
more  complicated  systems,  methods  development  is  an  important  component  of  our 
research  program.  To  this  end,  we  are  developing  new  techniques  for  treating  electron 
correlation  and  for  finding  transition  states  on  potential  energy  surfaces.  First,  we  have 
derived  equations  for  a  direct  full  configuration  interaction  algorithm  that  utilizes  the 
pseudospectral  approach.-4  To  date,  no  one  has  used  this  numerical  technique  for  Cl.  We 
have  shown  that  it  improves  the  scaling  of  full  Cl  from  CN4  to  CMN2,  where  C  is  the 
length  of  the  Cl  eigenvector,  N  is  the  number  of  molecular  orbitals,  and  M  is  the  number  of 
grid  points  used  in  the  pseudospectral  part  of  the  calculation.  Since  M  scales  like  N,  we 
achieved  a  factor  of  N  reduction  in  scaling,  which  will  allow  the  execution  of  full  Cl 
calculations  on  larger  molecules  than  has  been  possible  in  the  past.  We  have  verified  the 
accuracy  of  the  technique  by  comparing  the  total  energies  for  conventional  (spectral)  full 
Cl,  pseudospectral  full  Cl,  and  two  hybrid  spectral/pseudospectral  Cl  approaches.  We 
find  that  one  particular  hybrid  method  gives  total  energies  that  differ  on  the  order  of  1 
microhartree  from  the  spectral  results,  while  retaining  the  scaling  savings.  We  are  now 
moving  to  a  more  practical  objective,  namely  implementing  a  similar  algorithm  for  truncated 
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Cl  (e.g.,  multi-reference  single  and  double  excitation  Cl).  Second,  we  have  recently 
developed  an  entirely  new  approach  to  finding  saddle  points  on  potential  energy  surfaces. 
This  new  method  corresponds  to  reaction  barrier  following  rather  than  reaction  path 
following.  We  have  tested  this  on  both  empirical  and  ab  initio  potential  surfaces  and  find  it 
to  be  accurate,  elegant,  and  simple:  it  utilizes  the  energy  gradient  only  occasionally  and 
never  the  energy  second  derivatives  (except  at  the  very  last  to  check  the  rank  of  the  saddle 
point).  Further,  it  requires  no  initial  guess,  which  might  bias  the  solution,  for  the  saddle 
point  configuration.  Further  work  on  this  method  is  in  progress,  and  an  initial  paper  is  in 
preparation  describing  the  technique.25 
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Theoretical  Chemistry  Conference,  Salt  Lake  City,  Utah  (with 
Christine  Wu). 

"Study  of  the  Dynamical  Stucture  of  the  Reconstructed  Si(100) 
Surface  Via  Molecular  Dynamics  Simulation,"  talk  at  the  West  Coast 
Theoretical  Chemistry  Conference,  Salt  Lake  City,  Utah  (with 
Paul  Weakliem). 


"Fluorine  Atom  Etching  of  the  Silicon(  1  ()0)-2x  1  Reconstructed 
Surface,"  poster  at  the  American  Conference  on  Theoretical 
Chemistry.  San  Diego,  California  (with  Christine  Wu). 

"Study  of  the  Reconstruction  and  Structure  of  Semiconductor  (100) 
Surfaces,"  poster  at  the  American  Conference  on  Theoretical 
Chemistry,  San  Diego,  California  (with  Paul  Weakliem). 

"Structure  and  Etching  of  Group  IV  Semiconductors  by  Computer 
Simulation,"  talk  at  the  First  Annual  Surface/Interface  Research 
Meeting  of  the  Southern  California  AVS,  Riverside,  California  (with 
Paul  Weakliem). 

"Hydrogen  Adsorption,  Diffusion,  and  Desorption  on  the  Si(100)- 
2x1  Surface."  talk  at  the  West  Coast  Theoretical  Chemistry’ 
Conference,  Moffett  Field.  California  (with  Christine  Wu). 

"Efficient  Algorithms  for  Correlation-Consistent  Configuration 
Interaction,"  poster  at  the  West  Coast  Theoretical  Chemistry 
Conference,  Moffett  Field,  California  (with  Todd  Martinez  and 
Aseem  Mehta). 

"Theoretical  Studies  of  Hydrogen  Adsorption,  Diffusion,  and 
Desorption  on  the  Si(l(X))-2xl  Surface."  poster  at  the  Gordon 
Conference  on  Surface  Chemistry  of  Electronic  Materials,  Oxnard, 
California  (with  Christine  Wu). 
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"Initial  Stages  of  Fluorine  Etching  of  Si(100)  by  Molecular 
Dynamics  with  a  First  Principles  Many  Body  Potential,"  poster  at 
the  Gordon  Conference  on  Surface  Chemistry  of  Electronic 
Materials,  Oxnard,  California  (with  Paul  Weakliem). 

"Monte  Carlo  Simulations  of  SixGei-x/Si  and  SixGei-x/S=Ge 
Superlattice  Structures,”  poster  at  the  13th  Annual  West  Coast 
Theoretical  Chemistry  Conference,  Richland,  Washington  (with 
Lawrence  Carter). 

"H2  Desorption  and  H  Adatom  Diffusion  on  Si(100)-2xl,"  poster  at 
the  13th  Annual  West  Coast  Theoretical  Chemistry  Conference, 
Richland,  Washington  (with  Christine  Wu). 

"Pseudospectral  Full  Configuration  Interaction,"  talk  at  the  13th 
Annual  West  Coast  Theoretical  Chemistry  Conference,  Richland, 
Washington  (with  Todd  Martinez). 

"New  Ab  Initio  Methods  for  Searching  for  Transition  States,"  talk 
at  the  13th  Annual  West  Coast  Theoretical  Chemistry  Conference, 
Richland,  Washington  (with  Irina  Ionova). 
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